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Hamiltonian time evolution in terms of an explicit parameter time is derived for general relativity, 
even when the constraints are not satisfied, from the Arnowitt-Deser-Misner-Teitelboim-Ashtekar 
action in which the slicing density a(x,t) is freely specified while the lapse N = ag 1 ^ 2 is not. The 
constraint "algebra" becomes a well-posed evolution system for the constraints; this system is the 
twice-contracted Bianchi identity when Rij = 0. The Hamiltonian constraint is an initial value 
constraint which determines g 1 ^ 2 and hence N, given a. 
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A minor change (ljg| in the Arnowitt-Deser-Misner 
(ADM) action principle leads to striking conse- 

quences for the understanding of general relativity in 
Hamiltonian form. Recent work on hyperbolic formu- 
lations of general relativity indicates UpT] that 
the freely specifiable quantity that determines the slicing 
of spacetime is not the lapse N but the "slicing density" 
a(x,t) — Ng^ 1 / 2 (where g = det</y is the determinant 
of the spatial metric). Altering the action principle to 
take this into account leads to a number of key results: 
1) equations of motion that are equivalent to = 0, 
not Gij = 0, 2) a Hamiltonian vector field that gener- 
ates time evolution even when the constraints are not 
satisfied, 3) a constraint algebra that is a homogeneous 
symmetric hyperbolic system that dynamically preserves 
the constraints, and 4) a new understanding of the Dirac 
algebra. 

Before beginning, let us fix notation and note a few 
elementary results. We work on a manifold £ x R with 
a "foliation-adapted" cobasis for the metric 

ds 2 = -N 2 {dtf + g l3 {dx l + (3 i dt)(dx j + ftdt). (1) 

Here, N is the lapse function (a space scalar) and (3 l is 
the spatial shift vector. Quantities with overbars denote 
spatial quantities, in particular R the spatial curvature 
scalar obtained from g^ and Vj the corresponding spa- 
tial covariant derivative. We also introduce the extrinsic 
curvature K%j and its trace H = K k- The momen- 
tum conjugate to the metric is a density of weight one, 



lapse (the "slicing density") a = Ng^ 1 / 2 = a(x,t) is 
freely specifiable while N is not |lq |. (The slicing density 
a is also used prominently in JlJ and |^Jl^].) Indeed, if 
one computes do log a = f(x, t) from a given a{x, t), then 
one finds 
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1/2 d a = d N + N 2 H = Nf, 



(2) 



L2) with/(x,t) = 
Combining this 



the equation of harmonic time slicing [ 
<9ologa acting as a "gauge source" |18]. 
with the 3+1 equation for H from the trace of (|10j) below, 
one obtains a quasi-linear wave equation for N . Every 
foliation is described by such a wave equation for some 
value of a. This wave equation for N |p|,p|p^| played a 
vital role in the development of our first order symmetric 
hyperbolic "Einstein-Ricci" system |^,^0|-|l2]| and reflects 
the built-in causality which comes from working with a. 
We conclude that TV, which determines the proper time 
as NSt between slices t — t' and t = t' +St, is a dynamical 
variable (cf. [0), closely connected to g 1 ^ 2 . N can also 
be seen to be determined from a and the Hamiltonian 
constraint, written as the generalized and completed con- 
formal "Lichnerowicz equation" |ll]- Q ; the explicit form 
in is the suitable one for present purposes. From this 
perspective, the Hamiltonian constraint plays its familiar 
role as an initial value constraint which determines g 1 / 2 
given a complete set of freely specified data ^0|. The 
important insight is that this in turn determines N from 
a, so that the Hamiltonian constraint does not fix the 
time but does fix the proper time r's rate with respect 



7r y = g l/2 {Hg %1 - K %3 ). The natural time derivative for to t: dr/dt = ag 1/2 = N along the normal do- 



evolution, do, acts in the normal future direction to the 
spacelike slice and is denoted by an over-dot. It is given 



by do = dt — £ a, where £n is the Lie derivative along 
the shift [3. ' 

We have found in work on hyperbolic formulations 
of the equations of evolution of general relativity which 
have no unphysical characteristics p|-|l~5j] that we must 
in essence use the Choquet-Bruhat "algebraic gauge" ||] 
to restrict the ordinary lapse N . The weight-minus-one 



Motivated by these findings, we alter the undetermined 
multiplier in the ADM canonical action principle from N 
to a. Using a has the affect of altering the Hamiltonian 
density from TL to 



ft = g 1/2 H = n ij TTij - ivr 2 - gR, 



(3) 



the latter being of scalar weight plus two and a ratio- 
nal function of the metric. (Note we reserve the phrase 
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"Hamiltonian constraint" to refer to the equation H = 
and use "Hamiltonian density" for Tt, which may not 
vanish, similarly for the momentum constraint and den- 
sity.) This leads to Teitelboim's ||] and Ashtekar's ||] 
modification of the ADM action (IQttG = 1 = c) 



Sadmta[s, tt; 01,(5) = / d A x (n lj gij - aH) 



(4) 



where we use Kuchaf 's notation indicating functional and 
explicit function dependence and boundary terms are 
ignored here as they are not the focus of the present 
analysis. [There are no difficulties in obtaining, for ex- 
ample, the ADM energy surface integral. One requires 
the same asymptotic conditions as always, maintaining 
N -> 1 + 0{r~ x )— not a -> 1 + 0(r _1 )— while recall- 
ing N = g 1 / 2 a.} The vacuum case is considered, but 
to add minimally coupled matter and/or a cosmological 
constant is straightforward. The momentum density TLi 
has been absorbed into 7r lJ g^ through use of the time 
derivative 8q. Explicitly, the Lie derivative term in tt 4 - 7 
is, up to a divergence, 



(5) 



Consider a general variation of the modified Hamilto- 
nian density 

6H = (2nj - g l3 ir)5Ti VJ + (27T 4 V fc - 7nr ij (6) 
+gR ij - gg ij R)S 9ij - fl (V*V^ff« - g ij V k V k S 9ij ). 

Note that (^J) does not involve either the Hamiltonian 
or momentum densities while, in contrast, the variation 
of the ADM Hamiltonian density SH = Sig^^H) does 
contain a term proportional to the Hamiltonian density. 

Requiring that Sadmta be stationary under a varia- 
tion with respect to Tt lJ gives the definition of the extrin- 
sic curvature 



fJU 



sn 



a(27Ty - gijn) 



-2NKij. (7) 



Requiring that it be stationary under a variation with 
respect to </y gives the equation of motion 



5H 



Sgij 

+g(V i V j a-g ij V k V k a) 



ag{R i: > - g v R) - a(2TT ik ir j k - mr lJ ) 



(8) 



The slicing density a and the shift (3 Z are not to be 
varied. Instead the constraints arc imposed on initial 
data and are preserved dynamically as shown below. This 
is not an already parametrized theory in the usual sense. 

Consider the familiar 3+1 identities 



9ij = -2NK; 



131 



(9) 



K, 



V; />', , — Rij 



HK %3 - K lk K k - N-^iVjN). 

(10) 



Also, recall the formula for the derivative of the deter- 
minant of the 3-metric, g^g — g v - ' g-ij = —2NH. Now 
we pass to canonical variables. Using (||), ([To|), the time 
derivative of tt ij is computed to be identically 

7r y ' ee Ng 1 / 2 (Rg il - R lj ) - Ng-^ 2 ( 



+g 1/2 {V l V J N - g lJ V^V^TV) + Ng^'^TZ^} 



{/TT 7T j k 
1/2 r 



') 
(ID 



where IZij = Rij 



We see that the equations of motion (0) and (||) de- 
rived from the action principle are (^j and ( |TT| ) when 
i? u — g VJ R k k = 0. Thus, to say that (|) holds is to assert 
that R?i = 0. The equations of motion hold strongly, in- 
dependent of whether the constraints are satisfied. This 
is not true in the ADM formulation because of the pres- 
ence of the Hamiltonian density in their equation of mo- 
tion for 7!"^ . 

This difference can be explained more fully as follows. 
From the definition of the Einstein tensor in terms of the 
Ricci tensor, = R^ — ^g^R^a, and the observation 
that 2G°o = R°o — R k k i one derives the identity 



G, 



gijG°o = R^ — g%jR K 



(12) 



The vanishing of the right hand side does not depend 
on either the Hamiltonian or momentum densities and 
is equivalent to Rij — 0. Clearly it is also equivalent 
to Gij = —gijG°Q. Thus, while i? M „ = and G M „ = 
are equivalent, R^ — and Gij — are not equivalent 
as equations of motion — unless the Hamiltonian density 
Ti. = 2(? 1 / 2 G o vanishes exactly, that is unless the Hamil- 
tonian constraint holds. The ADM action principle is 
equivalent to Gij = and so also requires H — to be 



equivalent to ([Tlj). We recall that the use of Rij has al- 
ways been preferred by the French school, pioneered by 
Lichnerowicz |2^] and Choquet-Bruhat pi| . 

This raises an important principle: A constrained 
Hamiltonian theory should be well-behaved even when 
the constraints are violated. As discussed in |^|| , recent 
efforts fjL|l5| to achieve well-posed hyperbolic formula- 
tions of general relativity, with only physical character- 
istics, can be understood in this light as well. From this 
point of view, the Hamiltonian and momentum densi- 
ties are definite fixed combinations of the phase space 
variables, but their values may deviate from zero. The 
form of the equations of motion should not depend on 
these values. When the constraints are satisfied (i.e., 
the densities vanish), one is on the so-called constraint 
hypersurface, and many unphysical degrees of freedom 
are frozen because relations among many of the variables 
(g, 7r) are fixed. When the constraints are relaxed, the 
theory explores phase space away from the constraint hy- 
persurface. The objective is to have a theory whose char- 
acter does not change dramatically when one moves off 
the constraint hypersurface. Examples for which this is 
particularly relevant are numerical applications where vi- 
olation of the constraints is inevitable. It seems that sim- 
ilar properties may be shared by the canonical Ashtekar 
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variables, but there are subtleties beyond the scope of 
this paper that require closer investigation (cf. |l7| ). [At 
this point, we should stress that the Rij = equations 
given by (^) , (||) or @ , (|o|) are not in themselves known 
to be well-posed, though they have no unphysical charac- 
teristic speeds 0. They do, however, lead to well-posed 
evolution of H and Hi as we shall see.] 
Introduce the smeared Hamiltonian 



Ri 



as follows. Combine the identity (112) with the 



H a = / dV a(x' ,t)H. 



(13) 



The equation of motion for a general functional on phase 
space F[g, tv; x, t) is 

doF[g, x, t) = ~{H a , F[g, tt; x, t)} + d F{g, tt; x, t). 

(14) 

Here, 8q is a total time derivative while do is a "partial" 
derivative of the form 9 t — £ ^ which only acts on explicit 
spacetimc dependence. The Poisson bracket is given by 



{F,G} 



d A x 



5F 



SG 



SF 



SG 



Sgij {x, t) Stt 1 ^ (x, t) Stt^ (x, t) Sg^ (x, t) 



It is evident that the equations of motion (^), (@) are ob- 
tained by applying (|l4|) to the canonical variables (g, 7r). 
Also, observe that applying (jlj) to N — ag 1 / 2 produces 

>!• 

Time evolution is generated by the Hamiltonian vector 
field 

X H a = fd 3 x{a(2n lJ - ng l3 )^— - [ag(R l > - g ij R) (15) 

c 

+a(2n ik Tr j k - irn ij ) - g(V l V j a - g i3 V k S7 k aj\— — }. 

Sit 13 

Again, this does not depend on the Hamiltonian or mo- 
mentum densities, so it is a good time evolution operator 
even away from the constraint hypersurface. (This ob- 
servation was essentially made by Ashtekar in footnote 
17 of |Q but evolution was mistakenly associated with 
Ga = 0). 

By the product rule shared by do and the Poisson 
bracket, one computes the evolution equations for the 
constraints to be 

H = -{H a , H} = gag^diHj + 2gg i *H i V j a, (16) 
Hj = -{H Q , Hj) = ad 3 H + VHdjCL, (17) 

where V^a = dja + ag~ 1 ^ 2 dig 1 ^ 2 . These equations cor- 
respond to ( p0| ) and (^TJ) below when R^ — 0. Thus, the 
Poisson brackets of the smeared Hamiltonian with the 
unsmeared densities are seen to be well-posed evolution 
equations for the densities. 

These equations can be shown to be equivalent to the 
twice-contracted Bianchi identities V ' pG^ a = when 



twice-contracted Bianchi identities to obtain a transpar- 
ent form of the Bianchi identities 



VpG^o = V G°o + V ; G J o = 0, 



(18) 



V G^ = V G\ - V,-G o + Vil&j - 5 l R\] = 0. 

(19) 

(After the expression in square brackets is replaced using 

the vanishing of Rij, these become equations of motion 

rather than identities.) Express the Bianchi identities 

©, © in 3+1 language, using C = g^H = 2G° and 
r~t- „=TV2-i 



C 



T/2 Hi = 2NG°i. A calculation yields 



C - NV 3 Cj = 2 [CjV 3 N + NHC - NK lj [TZ i3 }j , (20) 
Cj - NVjC ee 2 (cVjN + \nHC 3 - V\N[TZ l3 })) , (21) 



where IZij = Rij 



ijR k- This system is clearly sym- 
metric hyperbolic with only the light cone as character- 
istic, and changing to H, Hi gives ( |l^) and (p^) when 
Rij = 0. (Without considering our identities, or a Hamil- 
tonian framework, Frittelli |26| reached the same conclu- 
sion about well-posedness of constraint propagation in 



the "standard" 3+1 formulation |19f|, which uses Ri 



0, 



and its absence in the ADM equations, with Gy = 
0. Related formulae were also obtained by Choquet- 
Bruhat and Noutchegueme [E7J for the evolution of mat- 
ter sources p 00 , p oi , where pP™ = T 13 " - \gP a T^ M .) 

Return to the Hamiltonian formulation. The Poisson 
bracket between two smeared Hamiltonians is 



{H ai ,H a2 } = -J d 3 xg{ ai V l a 2 - a^ax)^ (22) 

= — J d 3 x gg % 3 {ct\djd2 — a 2 dja\)Hi 

This bracket expresses the consistency of time evolution 
under different choices of a. The Jacobi identity is 

{H ai , {Ha 2 , F}} — {H a2 , {H ai ,F}} = 

= {{H ai ,H a2 },F}. (23) 

Because of the metric dependence in (H|), one sees that 
the difference between evolution with a 2 followed by a± 
and the reverse is a spatial diffeomorphism when Hi = 
@ (or when SF/Stt^ = 0). 

These results lead to a new understanding of the Dirac 
"algebra" of the constraints (cf. p8||). As is well-known, 
the Dirac algebra is not the spacetime diffeomorphism 
algebra. The root of this is that the action (^) is invari- 
ant under transformations generated by the constraints 
even when they are not satisfied p9[ . The equations 
which hold even when the constraints are not imposed 
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are R 



0. These equations are preserved by spa- 



tial diffeomorphisms and time translations along their 
flow, yet a general spacetime diffeomorphism applied to 
Rij = mixes in the constraints. A comparison of (16), 
@ and @, @ shows the effect clearly. The Bianchi 
identities are spacetime diffeomorphism invariant while 
the constraint evolution equations derived from the ac- 
tion principle are not. The equations ([l6|), jl7| ) and (20), 
( ^l| ) differ precisely by terms proportional to Rij . 

A second crucial understanding is the way in which 
the once-smeared form of the Dirac algebra ( |l6| ) , (JTt]) en- 
sures consistency of the constraints via a well-posed ini- 
tial value problem. If the constraints vanish initially, then 
they always vanish in a corresponding physical domain of 
dependence. This dynamical mechanism for consistency 
follows from the dual role of Ti. as part of the generator 
of time translations and as an initial value constraint. 

It is worth re-emphasizing the altered role of the 
Hamiltonian constraint. The Hamiltonian constraint is 
an initial value constraint from which g 1 / 2 is determined 
as in the solution of the initial value problem jljj , which 
then allows N to be reconstructed from a. By virtue of 
(plf), (fl7|), once the initial value problem is solved, it re- 
mains solved in a spacetime domain dictated by causality. 
The Hamiltonian constraint does not express the dynam- 
ics of the theory; ( p^ ) is the dynamical equation. The 
application of these ideas to canonical quantum gravity 
will appear elsewhere |}0| . 
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